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Abstract. We report the identification of the cataclysmic variable 1RXS J062518. 2+733433 as an intermediate 
polar. The orbital period of the system is determined to be 283.0 ± 0.2 min from the radial velocity variation of 
Ha, measured in an extensive set of time-resolved spectroscopy. Differential optical photometry obtained over a 
base line of three weeks reveals the presence of coherent variability with a period of 19.788 ± 0.002 min, which we 
suggest to be the white dwarf spin period. The power spectrum of our photometry also contains a strong signal 
near the spectroscopically determined orbital period. The emission lines in 1RXS J062518. 2+733433 display a 
complex multicomponent structure. In the trailed spectrogram of He I A 6678 we detect a narrow component with 
a radial velocity semi-amplitude of ~ 140 kms, consistent with a possible origin on the irradiated face of the 
secondary. The absence of eclipses gives an upper limit on the binary inclination of i 60°. 
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1. Introduction 

Over the past few years it has become increasingly clear 
that our understanding of the evolution of cataclysmic 
variables (CVs) is a rather unsettled matter. A number of 
the predictions made by the stand ard theory of CV evolu- 
tion (disrupted magnetic braking, lKinJl988j) are in direct 
conflict with the observations. Whereas a number of ad- 
ditions / alternati ves to the standard model have been re- 
cently proposed llKing fc Schenkerl2~0 02: Schcnke r fc Kmel 
l2002HSchenker et al.l2002HAndronov et al.l2003|) . it is ap- 
parent that the disturbing disagreements between the the- 
ory and the observations have a common denominator: the 
possible impact of selection effects on the currently known 
population of galactic CVs. In order to quantitatively test 
any theory of CV evolution it is mandatory to establish a 
statistically complete sample of this class of binary stars. 

There are currently a number of large scale surveys 
for CVs underway that p ursue this aim (e.g. lSzkodv et alJ 
120021: iMarsh et al.ll2002() . Among these projects, our se- 
lection of CVs based on their spectro scopic properties in 
the Hamburg Quasar Survey (HQS; lHagen et alJ IT995) 
has been especially prolific, resulting in the discovery of 
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more than 50 new bright CVs fe.g iGansicke et alJl200ol 
iNogami et aT1l200ll iGansicke et al.ll2002alblh ~ 

In this paper, we report follow-up observations of 
1RXS J062518.2+733433, henceforth RX J0625, an object 
originally identified as a C V on the base of i ts X-ray emis- 
sion and optical spectrum l|Wei et alJllQ QQ!). We have in- 
dependently selected RX J0625 as a CV candidate because 
of the noticeable Balmer emission in its HQS data. The 
strong He n A 4686 emission detected in the identification 
spectrum of RX J0625, along with coherent optical vari- 
ability observed on a time scale of ~ 20 min, immediately 
reveal the intermediate polar nature of this CV. In Sect. 2 
we describe the observational data and its reduction. We 
then derive the orbital period from the radial velocity vari- 
ation of the emission lines and the spin period from differ- 
ential photometry in Sect. 3. Finally, in Sect. 4 we discuss 
the behaviour of the strongest emission lines in detail and 
summarise our findings in Sect. 5. 



2. Observations and Data Reduction 

2.1. Spectroscopy 

On 2001 April 29 we obtained a single identification spec- 
trum of RX J0625 at the Calar Alto 2.2m telescope with 
the CAFOS focal reductor spectrograph, using the stan- 
dard SITe CCD (Table 1). We used the B-200 grism and a 
slit width of 2", resulting in a useful wavelength range of 
3500-7000 A and a spectral resolution of 4.7 A (Fig.P). A 
flux standard (BD+75 325) was observed with the same 
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Table 1. Log of Observations. 
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Fig. 1. Identification spectrum of RX J0625 obtained on 
2001 April 29. 



set up in order to correct for the instrumental response. 
The identification spectrum contains noticeable Balmer, 
He i and He u and N/C Bowen emission lines. The strength 
of He ii A 4686 is comparable to H/3, indicating the pres- 
ence of a strong source of ionising photons in RX J0625, 
typical of either magnetic CVs or novalike variables. In 
order to determine the orbital period of RX J0625, we ob- 
tained 45 higher resolution spectra, again with CAFOS 
at the Calar Alto 2.2m telescope (Table 1). This time the 
G-100 grism was used in conjuction with a 1.2" slit which 
gave a wavelength range of 4240-8300 A and a spectral 
resolution of 2.1 A. The higher resolution spectra were ob- 
tained over a period of 3 weeks, optimising the sampling 
for an efficient period determination. 

All spectra were reduced in a standard manner using 
the Figaro package within the Starlink software collec- 
tion. The frames were corrected for the bias level by sub- 
tracting the mean of a series of bias images taken at the 
start and end of each observing night. Dome flat-fields 
were used to remove pixel to pixel variation s of the chip . 
The spectra were then optimally extracted (Home 1986) 
and sky line subtracted using Tom Marsh's Pamela pack- 
age. Especial care was taken to account for the tilt of 
the spectra in order to maximise the signal-to-noise ratio. 
The wavelength calibration was performed using mercury- 
cadmium, helium-argon and rubidium arcs. Uncertainties 
on every data point calculated from photon statistics are 
rigorously propagated through every stage of the data re- 
duction. We did not attempt to flux-calibrate the higher 
resolution spectra as it was not required for the radial 
velocity analysis described below. 

Table 2 gives the equivalent widths and line widths of 
the main emission lines detected in the average of the 45 
high resolution spectra. 
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Table 2. Equivalent widths and line widths (corrected for 
the instrumental resolution) of the strongest emission lines 
measured from the average of the 45 spectra obtained in 
December 2002. 
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2.2. Photometry 

We have obtained differential U-band and filter-less pho- 
tometry of RXJ0625 during 6 nights in December 2002 
with the CAFOS SITe CCD camera on the Calar Alto 
2.2m telescope (TableQ]). In order to achieve a high time 
resolution (~ 30s), only a small window of the chip 
(~ 3' x 2') was read out. The data were bias-subtracted 
and flat-fielded in a standard fashion using the ESD-MIDAS 
package, an d aperture magnitudes w ere extracted with the 
sextractor l|Bertin fc Arnoutslll996|) . The V magnitudes 
of RX J0625 were derived relative to the HST Guide Star 
GSC 0437000234 (V = 13.4 ± 0.4), located 1'50" south- 
wards of RX J0625. This comparison star has been satu- 
rated in a number of the filter-less CCD images obtained 
on December 14/15, and, hence, we used GSC 0437000998 
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Fig. 2. Differential CCD F-band (V) and filter-less (C) 
photometry obtained at the Calar Alto 2.2m observatory. 
Note that the different scales of the filter-less and V band 
data. 
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(V = 14.9±0.4), located l'lO" northwards of RXJ0625 for 
the reduction of these images. RX J0625 was found at an 
average magnitude of 14.80 ± 0.05 throughout the nights 
in which we used the V band filter. The light curves of 
RXJ0625 clearly reveal the presence of variability with 
a period of ~ 20min and an amplitude of ~ ±0.2 mag 
throughout all nights (FigJ5J. In addition to this short- 
term variations the light curve of RXJ0625 displays a 
modulation on time scales of several hours. 



3. Analysis 

3.1. Spectroscopy 

The primary aim of our time resolved spectroscopy of 
RX J0625 is to measure its orbital period from the radial 
velocity variations of the emission lines. The line profiles 
clearly display a complicated multicomponent structure 
(Fig-El - I n a fi rs t attempt, we used a single Gaussian least 
square fitting procedure to determine the radial velocity 
variations of Ha, H/3 and Hen A 4686. The radial velocity 
data were then subjected to the following pe riod analysis 
meth ods: (a) the phase dispersio n method llStellingwerJ 



Fig. 3. Continuum-subtracted and normalised Ha profile 
sample sorted by orbital phases computed from Eq. . 



1978 ) , (b) the analysis of variance ijSchwarzenberg-Czernv 
1989, as implemented in the MIDAS context TSA), (c) and a 



sine wave x 2 fitting procedure, using chisq in the Starlink 
package period. The resulting periodograms for Ha are 
shown in Fig.|3] The orbital frequency (period) is inferred 
from the three independent techniques in a consistent way 
to be Q = 5.088±0.003d- 1 (P orb = 283.0±0.2min), where 
the error is computed from the \ 2 sme wave fitting (2a, 



108 6 4 
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Fig. 4. Period analysis of the radial velocities measured 
from a Gauss-fit to the Ha emission line profiles, (a): 
Phase dispersion method periodogram. (b): Analysis of 
variance periodogram. (c): % 2 sine fit. 
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Fig. 5. Radial velocities measured from a Gauss- fit to the 
Ha emission line profiles, folded with the orbital period 
Porb = 283.0 min, using Eq. (J2J). The dashed curve is the 
best sine fit to the radial velocities. 

4 interesting parameter). A slightly more conservative er- 
ror estimate based on the results obtained from all three 
methods gives P or b = 283.0 ± 0.3 min. The analysis of H/3 
and He ii A 4686 provides consistent results for Porb, but 
with significantly larger errors. The robustness of our pe- 
riod determination is shown in Fig. [S] were we folded the 
Ha radial velocities over the orbital period. The radial ve- 
locity measurements from this fit are provided in Table 3. 

For comparison purposes we measured the radial ve- 
locities of Ha cross-correlating the observed line pro- 
file with a single Gaussia n of fixed width of 300 km s -1 
( Sch neider &: Youndll980j) . The resulting orbital period is 
identical to that derived from the single Gaussian fit to 
the Ha profiles, but the radial velocity variation displays 
an amplitude larger by a factor of ~ 2. 

Finally, we attempted to fit the Ha emission with a 
blend of a narrow and a broad Gaussian component to 
model more accurately the complex structure of the ob- 
served profiles. Unfortunately, this procedure did not pro- 
vide unambiguous results because of the limited spectral 
resolution of our data. Due to the multicomponent struc- 
ture of the lines profile (see Fig. 02) it is difficult to in- 
terpret the radial velocity variation obtained from either 
the single Gaussian fit (FigJSJl or the Gauss correlation to 
Ha. The behaviour of the most important emission lines 
is discussed in more detail in Section the basis of trailed 
spectrograms. 

3.2. Photometry 

In order to analyse the periodicities observed in the light 
curve of RXJ0625 (Fig.|2Jl, we have computed fro m the 
entire photometric data set both a IScargld fl982) peri- 
odogram (as implemented in the MIDAS context TSA) as 
well as a p eriodogram using t he Phase Dispersion Method 
(PDM) of lStellingw erf (1978). Because the lengths of our 
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Fig. 7. Period analysis near the presumed white dwarf 
spin. Top panel: Scargle periodogram. Bottom panel: 
phase dispersion periodogram. The black lines are com- 
puted from the observed data. The gray lines are com- 
puted from a sinusoid assuming a period of 1/w = 
19.788 min and identical sampling as in the observed data. 



individual photometric observations are of the order of 
or shorter than the orbital period subtracting the nightly 
mean from the data would introduce erroneous signals in 
the low-frequency range of the periodograms. Considering 
that RXJ0625 does not exhibit noticeable night-to-night 
variability, we therefore decided to subtract from the V 
filter and white light data the mean of all measurements 
obtained in the corresponding band prior to the compu- 
tation of the periodograms. The resulting periodograms 
(Fig-EI an d d contain strong signals concentrated in the 
frequency ranges / Si 10 d -1 and / ~ 60 — 70 d -1 . The 
nature of the signals detected in these two separate fre- 
quency regimes is discussed below. 

High-frequency signal. Figure|7| shows an enlargement 
around the peak signal contained in the high-frequency 
range of the periodogram shown in Fig.El Both pe- 
riod analysis methods consistently detect the maximum 
power at 72.772 ± 0.008 d _1 , corresponding to a period of 
19.788 ± 0.002 min, whose error is given by la width of 
the Gaussian fit to the frequency peak. In order to test 
the significance of this signal, we have created a faked 
data set from a sine wave with a period of 19.788 min, a 
similar amplitude as the observed short-period variabil- 
ity, and sampled exactly in the same way as the observed 
data. In addition, we added a 10% error to these synthetic 
data. The periodograms for this faked data set are plot- 
ted in gray in Fig.0 and reproduce very well the entire 
structure of aliases contained in the periodograms of the 
observed data. We conclude from this comparison that 
the frequency of the maximum signal detected in the peri- 
odograms is indeed the true period underlying the short- 
term variability observed in RX J0625. 
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Table 3. Radial velocities of the Ha emission line determined from a single-Gaussian fit to the line profile (corrected 
for the heliocentric velocity). 
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Fig. 6. Scargle periodogram of the differential CCD photometry shown in Fig.0 u> and are the white dwarf spin 
and orbital frequency, respectively. 



Coherent variability on time scales of a several minutes 
to several tens of minutes has been detected in the opti- 
cal light curves of a number of intermediate polars (e.g. 
|Pattersonlll99^ . and is interpreted as the spin period of 
the accreting magnetic white dwarf. We conclude from the 
spectral appearance of RXJ0625 and from the detected 
coherent optical variability that RX J0625 is indeed a new 
member of the small class of intermediate polars. We sug- 
gest that the detected frequency (period) of 72.772 d -1 
(19.788 min) is the white dwarf spin frequency uj (period 
-Pspin = From the combined V band and filter- less 

photometry, we derive the following spin ephemeris: 

0o pin = HJD 2452617.0645(4) + 0.013742(3) x E (1) 

where spin — is defined as the spin pulse maximum. 
Errors in the last digit are given in brackets. Figurc|HIa) 
shows our photometric data folded with this ephemeris. 



The quasi-sinusoidal shape of this light curve, which is 
typical for the spin light curve of interm ediate polars (see, 
e.g., the spin light curve of FO Aqr by Ide Martino et all 
1994), lends support to our hypothesis that the 72.772 d 
frequency detected in the power spectrum of RX J0625 is 
indeed the spin frequency of the white dwarf. 

Another hallmark of intermediate polars is the detec- 
tion of beat frequencies between the white dwarf spin fre- 
quency u> and the orbital frequency fi, which arise from 
the reprocession of X-rays emitted from close to white 
dwarf surface on, e.g., the secondary star. Such sideband 
signals have been detected i n different sys tems atw- 20, 
uj — f2, uj + f2, uj + 2f2 (see IWarnerl|l986l for an interpre- 
tation of these frequencies). In RXJ0625, the only un- 
ambiguous detection of a sideband signal is uj — f2 (see 
Fig. E| . The photometric data folded over this beat period 
of 21.275 min are shown in Fig.|HJb). While we are confi- 
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dent that our identifications of the spin and beat period 
are correct we mention as a note of caution that time- 
resolved X-ray and/or polarimetric data are necessary to 
unambigously confirm this interpretation. 

Low-frequency signal. The strongest signal in the low- 
frequency range does not coincide with the orbital fre- 
quency derived in Sect . 13 .ll from the radial velocity varia- 
tion of Ha, but is found at 4.944 d _1 , corresponding to a 
period of 291.3 min. This period is ~ 3% longer than the 
orbital period. Considering the sampling of our photomet- 
ric time series, especially the fact that none of our data sets 
covers significantly more than one binary orbit, the most 
likely hypothesis is that the detected low-frequency sig- 
nal is related to the orbital motion of the binary. Further 
observations would be useful to test in detail whether the 
low-frequency photometric frequency is indeed identical to 
the orbital frequency O. 
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4. Discussion 

As already mentioned in Sect. 13. ll the line profiles clearly 
display a multicomponent structure (Fig.EJ, which makes 
the interpretation of the determined radial velocity varia- 
tions rather ambiguous (Fig.|5J| . To explore the behaviour 
of the emission lines in more detail, we constructed trailed 
spectra of Ha, H/3, He i A 6678, and Hen A 4686, which 
are shown in FigureESJ These diagrams were computed 
from the continuum-normalised spectra after binning into 
15 phase intervals. The He i A 6678 line clearly shows two 
emission components: a narrow one with a radial velocity 
semi-amplitude of « 140 km s" 1 and a wider one with ve- 
locity reaching ~ ±500 km s -1 . Both components are not 
in phase, the phase offset being ~ 0.15. The trailed spec- 
tra shown in Fig.Elwere phased according to the following 
orbital ephemeris 

0g rb = HJD 2452617.5083(9) + 0.1965(1) x E (2) 

where the instant of zero phase was derived as the 
time of blue-to-red crossing of the narrow component of 
He i A 6678. 

In analogy to the narrow emission lines observed in 
many magnetic CVs, this narrow emission line compo- 
nent may have its origin on the irradiated face of the 
secondary star. In order to test this hypothesis, we per- 
formed the following simple calculation. The secondary 
mass was estimated from t he m ass-period relation de- 
rived by ISmith fc DnilioTI jl998l) . which gives a value 
of M2 « 0.48 M Q . For the primary mass we adopted 
Mi = 0.8 Mq (the avera ge white dwarf mass in long- 
period CVs determined bv lSmith k Dhillonlll998|) . result- 
ing in a mass ratio q = M2/M1 ~ 0.6. Assuming a disc 
radius of 0.8i?Li, the lack of eclipses in the optical light 
curve of RX J0625 sets an upper limit on the orbital incli- 
nation of i ^ 60°. The projected velocity of the secondary 
is given by: 

= 2nasim 

Porb(l+<?) V ' 



Fig. 8. The photometric data points from Fig-EI 
mean-subtracted and folded over: (a) The spin period 
of 19.788min, using Eq. JTJ. (b) The beat period of 
21.275 min. (c) The orbital period of 283.0 min, using 
Eq. ©. 

where a is the binary separation. With the above assump- 
tions we the radial velocity variation of the secondary star 
can be as large as K 2 ~ 130 km s -1 , which is comparable 
to the value obtained from the trailed spectrogram of the 
narrow He 1 A 6678 line. The predicted K\ velocity is then 
« 80kms _1 , somewhat larger than the semi-amplitude of 
the radial velocity curve of the broad component of Ha 
(Fig-El see also Fig-EI, but not inconsistent. 

Taken at face value, the behaviour of the narrow emis- 
sion line components detected in our (admittedly rather 
low resolution) phase-resolved spectroscopy are consistent 
with an origin on the irradiated face of the secondary. 

As we mentioned above, the broad component of 
the He 1 A 6678 line displays a phase offset of ~ 0.15 
with respect to the narrow one, reaching a velocity of 
~ — 500 km s -1 at its maximum excursion to the blue. 
The phasing with respect to the narrow component (and 
assuming that our above interpretation of the origin of 
the narrow component is correct) indicates that the broad 
emission is not coming from an axisymmetric structure 
around the white dwarf. Nevertheless, the high veloc- 
ity dispersion observed in the broad component clearly 
points towards an origin close to the white dwarf. The 
maximum velocity to the blue is reached at phase ~ 
0.4 — 0.5. A similar high-velocity S-wave has been ob- 
serve d in the interme diate polars EXH ya ijHellier et alJ 
Il989|) and V1025Cen (jHellier et alJl2002|) . The same com- 
ponents are also characteristic of the SW Sextantis stars, 
which have been recently propo sed to be magnetic sys- 
tems l|R.odrfguez-Gil et al. 2001). In all these systems, it 
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Fig. 9. Trailed spectrograms of Ha, H/3, He i A 6678, He n A 4686, computed from the 45 spectra obtained in December 
2002 (Table 1). Orbital phases have been computed using Eq. 



is believed that these high-velocity S-waves form in the 
vicinity of the primary's magnetosphere. 

The trailed spectra of He n A 4686 seem to exhibit the 
same two emission components, with the narrow one also 
dominating and same velocity amplitudes as in He i A 6678. 
Also the trailed spectra of Ha and H/3 clearly show a mul- 
ticomponent structure, even though their narrow compo- 
nents are less obvious than in He i A 6678. In order to se- 
curely identify the different line components with physical 
emission sites in the binary, it will be necessary to obtain 
time-resolved spectroscopy study spanning several consec- 
utive binary orbits with better spectral resolution. 

5. Conclusions 

We have observed RX J0625 as part of our ongoing search 
for new CVs selected on the base of their spectroscopic 
properties from the Hamburg Quasar Survey. From the 
radial velocity variations measured in Ha, we determine 
an orbital period of 283.0 ± 0.2 min. The detection of co- 
herent optical variations clearly classifies RX J0625 as be- 
ing a member of the small class of intermediate polars. 
The period of these variations, 19.788 ± 0.002 min, which 
is most likely the spin period of the white dwarf. In terms 
of its orbital and spin period, RX J0625 is very similar to 
FO Aqr (P rb = 291 min, and P sp i n = 21 min). Our phase- 
resolved spectroscopy clearly shows that the emission lines 
are multicomponent structured, and we identify a nar- 
row component with a radial velocity semi-amplitude of 
~ 140 kms, which might originate on the irradiated face of 
the secondary. Extensive photometric monitoring, phase- 
resolved high-resolution spectroscopy and pointed X-ray 
observations of RX J0625 are strongly encouraged. 
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Note added in proof: After submission of this paper, we be- 
ca me aware of an indep endent photometric study of RX J0625 
by IStaude et alJ (I2003T) that confirms the likely white dwarf 
spin period and shows that the low-frequency photometric 
signal detected in RXJ0625 is consistent with our spectro- 
scopic orbital period within the errors of the measurements. 
Furthermore, the referee drew our attention to the fact that an 
orbital period for RX J0625 has been included on De cember 13 
2002 in the online catalogue of iDownes et alJ J2001), which is 
consistent with our value. 



References 

Andronov, N., Pinsonneault, M., & Sills, A. 2003, ApJ, 582, 
358 

Bertin, E. & Arnouts, S. 1996, A&AS, 117, 393 

de Martino, D., Buckley, D. A. H., Mouchet, M., & Mukai, K. 

1994, A&A, 284, 125 
Downes, R. A., Webbink, R. F., Shara, M. M., et al. 2001, 

PASP, 113, 764 

Gansicke, B. T., Fried, R. E., Hagen, H.-J., et al. 2000, A&A, 
356, L79 

Gansicke, B. T., Hagen, H. J., & Engels, D. 2002a, On the 
space density of cataclysmic variables, ed. B. T. Gansicke, 
K. Beuermann, & K. Reinsch (ASP Conf. Ser. 261), 190-199 

Gansicke, B. T., Hagen, H.-J., Kube, J., et al. 2002b, 
HS 0455+8315: A new eclipsing novalike variable, ed. B. T. 
Gansicke, K. Beuermann, & K. Reinsch (ASP Conf. Ser. 
261), 623-624 

Hagen, H. J., Groote, D., Engels, D., & Reimers, D. 1995, 
A&AS, 111, 195 



8 



S. Araujo-Betancor et al.: 1RXS J062518. 2+733433: A new intermediate polar 



Hellier, C, Mason, K. O., Smale, A. P., et al. 1989, MNRAS, 
238, 1107 

Hellier, C, Wynn, G. A., & Buckley, D. A. H. 2002, MNRAS, 
333, 84 

Home, K. 1986, PASP, 98, 609 
King, A. R. 1988, QJRAS, 29, 1 

King, A. R. & Schenker, K. 2002, in The Physics of 
Cataclysmic Variables and Related Objects, ed. B. T. 
Gansicke, K. Beuermann, & K. Reinsch (ASP Conf. Ser. 
261), 233-241 

Marsh, T. R., Morales-Rueda, L., Steeghs, D., et al. 2002, in 
The Physics of Cataclysmic Variables and Related Objects, 
ed. B. T. Gansicke, K. Beuermann, & K. Reinsch (ASP Conf. 
Ser. 261), 200-207 

Nogami, D., Kato, T., Baba, H., et al. 2001, MNRAS, 322, 79 

Patterson, J. 1994, PASP, 106, 209 

Rodriguez-Gil, P., Casares, J., Martmez-Pais, I. G., Hakala, P., 

& Steeghs, D. 2001, ApJ Lett., 548, L49 
Scargle, J. D. 1982, ApJ, 263, 835 

Schenker, K. & King, A. R. 2002, in The Physics of 
Cataclysmic Variables and Related Objects, ed. B. T. 
Gansicke, K. Beuermann, & K. Reinsch (ASP Conf. Ser. 
261), 242-251 

Schenker, K., King, A. R., Kolb, U., Wynn, G. A., & Zhang, 

Z. 2002, MNRAS, 337, 1105 
Schneider, D. P. & Young, P. 1980, ApJ, 238, 946 
Schwarzenberg-Czerny, A. 1989, MNRAS, 241, 153 
Smith, D. A. & Dhillon, V. S. 1998, MNRAS, 301, 767 
Staude, A., Schwope, A. D., Krumpe, M., Hambaryan, V., & 

Schwarz, R. 2003, A&A, submitted 
Stellingwerf, R. F. 1978, ApJ, 224, 953 

Szkody, P., Anderson, S. F., Agiieros, M., et al. 2002, AJ, 123, 
430 

Warner, B. 1986, MNRAS, 219, 347 

Wei, J. Y., Xu, D. W, Dong, X. Y, & Hu, J. Y. 1999, A&AS, 
139, 575 



